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Abstract 

We employ an improved estimator to calculate the constraint effec- 
tive potential of the staggered magnetization in the spin | quantum 
Heisenberg model using a loop-cluster algorithm. The first and second 
moment of the probability distribution of the staggered magnetization 
are in excellent agreement with the predictions of the systematic low- 
energy magnon effective field theory. We also compare the Monte Carlo 
data with the universal shape of the constraint effective potential of the 
staggered magnetization and study its approach to the convex effective 
potential in the infinite volume limit. In this way the higher-order low- 
energy parameter is determined from a fit to the numerical data. 
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1 Introduction 



Systematic low-energy effective field theory is a very powerful tool for investigating 
the low-energy physics of Goldstone bosons. In particle physics this tool has been 
used to describe the low-energy physics of pions [ij. Pions arise as the pseudo- 
Nambu-Goldstone bosons of the spontaneously broken SU{2)l x SU{2)r chiral 
symmetry of QCD. Chiral perturbation theory provides a systematic momentum 
expansion of the low-energy physics of pions Based on symmetry considera- 
tions, the effective theory makes detailed quantitative predictions, which depend 
on the values of some a priori unknown low-energy parameters, such as the pion 
decay constant and the chiral order parameter. Systematic effective field theories 
have also been developed for quantum field theories and classical spin systems with 
a spontaneously broken 0{N) symmetry Effective field theories with an 0(3) 
symmetry have been used to describe the low-energy physics of magnons — the 
Goldstone bosons of quantum antiferromagnetism jJlEl, @, 0]- In this case, the rel- 
evant low-energy parameters are the spin stiffness p^, the spinwave velocity c, and 
the staggered magnetization density M.s. In particular, the effective theory makes 
detailed predictions for the finite-size and finite-temperature effects of systems in 
large volumes at low temperatures. These effects have been worked out by Hasen- 
fratz and Niedermayer even at the 2-loop level [8|. By comparing the effective field 
theory predictions with Monte Carlo data obtained with an efficient loop-cluster 
algorithm [9|, the low-energy parameters p^, c, and AAg have been determined with 
high accuracy for the spin ^ quantum Heisenberg model, both on the square 



and on the honeycomb lattice |12|. as well as for the J-Q model with 2- and 4-spin 
interactions on the square lattice [l3|]. Similar methods have been apphed to the 4-d 
0(4)-symmetric quantum field theory describing the Higgs sector of the standard 
model By comparing with effective field theory predictions for the suscepti- 

bility of the order parameter, both the vacuum expectation value of the Higgs field 
and the Higgs boson mass have been determined from Monte Carlo data for the 



4-d lattice 0(4) model [15|, llGj] obtained with the Wolff cluster algorithm [17|. The 
systematic low-energy effective field theory also predicts the shape of the constraint 
effective potential of the order parameter, which has been worked out in great detail 



by Gockeler and Leutwyler |l8l. Il9[|. Their predictions for the shape of the constraint 



effective potential have been teste d ag ainst Monte Carlo simulations of classical 3-d 



0(3) and 4-d 0(4) lattice models [20 



In this paper, we develop an improved estimator for the probability distribution 
of the staggered magnetization in order to extract the constraint effective potential 
from loop-cluster simulations of the spin | quantum Heisenberg model on the square 
lattice. The first moments of the probability distribution are in excellent agreement 
with the predictions of the magnon effective theory. The Monte Carlo data approach 
the universal shape of the effecti ve p otential in the infinite volume limit, which is 



known to be a convex function [2ll]. Our study confirms that the effective field 



theory provides quantitative predictions for a wide class of low-energy observables. 
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which are exact, order by order in a systematic momentum expansion. 

The rest of this paper is organized as follows. Section 2 summarizes the predic- 
tions of the low-energy magnon effective field theory that are relevant for our study, 
and section 3 discusses the determination of the low-energy parameters. In section 
4 we describe the improved estimator that allows us to obtain very accurate Monte 
Carlo data for the probability distribution of the staggered magnetization. In sec- 
tion 5 the results of numerical simulations obtained with a loop-cluster algorithm 
are compared with the effective field theory predictions. Finally, section 6 contains 
our conclusions. 



2 Effective Field Theory Predictions 



In this section we summarize the predictions of Gockeler and Leutwyler [18|, 119 
which are derived from a 3-d 0(3)-symmetric effective field theory. Although these 
authors had applications to the 3-d classical Heisenberg model in mind, their results 
also apply to the antiferromagnetic (2 + l)-d quantum Heisenberg model. This 
microscopic model is defined by the Hamiltonian 

H = jJ2s.-Sy-Ms-Bs, (2.1) 

where x and y are nearest-neighbor sites on a square lattice with spacing a. The 
spin i operators Sx obey the standard commutation relations 

[Sx, Sy] = iSxySabcSx- (2.2) 

— * 

Note that we work in natural units in which h = 1. Furthermore, Bg represents an 
external staggered magnetic field and 

Ms = ^(-l)(^i+^2)/'^S; (2.3) 

X 

is the staggered magnetization order parameter. In the infinite volume limit at zero 
temperature and with Bg = 0, Mg develops a non-zero vacuum expectation value 
signaling the spontaneous breakdown of the SU (2)^ spin symmetry down to its U{l)s 
subgroup. 

The corresponding low-energy effective field theory is formulated in terms of the 
staggered magnetization order parameter field 

e{x) = iei{x), e2{x), e^ix)) E S\ e{x f = 1, (2.4) 

where x = {xi,X2,t) is a point in Euclidean space-time. The leading terms in the 
effective action for the staggered magnetization field take the form 



S[e\ = J d?x dt 



^ (die- die + \dte- dtc] - M^e- 5, 



(2.5) 
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where ps is the spin stiffness, c is the spinwave velocity, and M.s is the staggered 
magnetization density. The corresponding partition function is given by 

Z = [ De exp{-S[e\). (2.6) 



Antiferromagnetic magnons have a "relativistic" dispersion relation with the spin- 
wave velocity c playing the role of the velocity of light. In fact, by introducing 
X3 = ct, the effective action can be written in the manifestly Euclidean rotation- 
invariant form 

S[^ = jd'x^ (^d^e-d^e-Mse-Bs^ , (2.7) 

which ensures Lorentz-invariance after analytic continuation from Euclidean to real 
time. It should be pointed out, however, that Euclidean rotation invariance is 
just an accidental symmetry of the leading terms of the effective action. Since the 
underlying quantum Heisenberg model itself does not have this symmetry, some 
of the higher-order four-derivative terms in the effective action break Euclidean 
rotation-invariance down to the discrete rotation subgroup of the square lattice. 

It is interesting to note that the ferromagnetic quantum Heisenberg model (which 
differs from the antiferromagnet only by the sign of the Hamiltonian) has very differ- 
ent symmetry properties at low energies. Unlike quantum antiferromagnets, quan- 
tum ferromagnets have a conserved order parameter — the uniform magnetization. 
Consequently, ferromagnetic magnons have a nonrelativistic dispersion relation and 



the corresponding effective action contains an additional Wess-Zumino term [22 
which breaks rotation invariance between space and Euclidean time already at lead- 
ing order. The resulting effective field theory for ferromagnetic magnons has been 



studied in detail in 23 



Here we concentrate entirely on antiferromagnets which are described by the 
effective action of eq. fl2.5p . We consider the system in a periodic cubic space-time 
volume L X L X (3 with the inverse temperature fixed aX (3 = L/c and with Bg = 0. 
The space-time average of the staggered magnetization is given by 



In contrast to [18|, ll9|] , we have included a factor | in the definition of $ because the 
quantum spins of the underlying Heisenberg model have S = \ while the effective 
field e{x) is normalized to 1. Due to the 0(3) symmetry, the probability distribution 

vm = \j^^ exp(-5[e1) 5[$-\^^j d'x e(x)^ (2.9) 

of the mean staggered magnetization vector $ only depends on the magnitude $ = 
|$|. The distribution is normalized by 

POO 

An / d$ $2 p($) = 1. (2.10) 
Jo 
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The constraint effective potential represents the free energy density of config- 
urations constrained to a fixed mean staggered magnetization $ and is determined 
by 

p{<^) =Arexp{-L\{<^)), (2.11) 



f2.12) 



with the normahzation derived in given by 

1 psL 



AT 



Svr^ce^o 



Here Aig = A^sC^^ is the staggered magnetization per spin, while (3q and /?i are 
shape-coefficients of the space-time box. For the exactly cubic space-time volume 
considered here jSo = 1.45385 and jSi = 0.225785. In the infinite volume and zero 
temperature limit the constraint effective potential approaches the infinite volume 



effective potential which is known to be a convex function of $ 211]. In a finite 



volume, on the other hand, m($) is not convex. An extensive variant of the intensive 
quantity «($) is 

U{^)=L^u{^). (2.13) 



Gockeler and Leutwyler have used chiral perturbation theory to systematically 
work out the finite-size effects of the constraint effective potential near its minimum 
(I8|,ll9| 



f/($) = UoW + 



(2.14) 



Here Uo{ip) and Ui{iIj) are L-independent functions of the rescaled variable 



(2.15) 



The leading order contribution to the constraint effective potential is given by an 
inverse Laplace transform 



dx exp{—ixip + T{ix)) 



of the function 



r(x) = Yl 



fSfiX 



?i=0 



(2.16) 



(2.17) 



Here the f3n are shape-coefficients of the finite space-time box described in detail 
in appendix B of (3|. As was pointed out in 19|, the function T{ix) is entirely 
kinematical and thus, unlike f/i(^/'), the quantity Uq{iP) is universal, i.e. completely 
independent of the low-energy parameters. Consequently, Uo{ip) is the same for all 
3-d systems with an 0(3) symmetry spontaneously broken down to 0(2), including 
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the 3-d classical and the (2 + l)-d quantum Heisenberg model. The 1/L correction 
to the leading contribution t/o('0) is given by 

Ui{'4)) = ^jj + ko exp{Uo{ijj)) / dx x"^ Re[exp{—ixip + T{ix))] 



#2 



dtp 



{2.18) 



Here is a low-energy constant related to the higher-order terms 



AS[e\ 



d X dt 



in the effective action, which is given by 



kn 



{hi + h^). 



(2.19) 



(2.20) 



It should be noted that the value of fco has no impact on eq. fl2.12p because it affects 
the normalization factor M only at higher orders in 1/L. 

Just as a non-zero quark mass in QCD explicitly breaks chiral symmetry and 
gives the pion its mass, a non-zero staggered magnetic field Bs explicitly breaks 
the SU{2)s spin symmetry and turns the magnons into pseudo-Nambu-Goldstone 
bosons with a non-zero mass m determined at leading order by 



MsBs 

PsC^ 



5. 



\BA. 



(2.21) 



The constant fco also appears in the 5^- dependence of the field expectation value 



c 
psL 



2^— n-("^cL) 



,n=0 



{mcLy 



+ ko 



mc 



Ps 



+ 0{nf 



(2.22) 

It should be noted that eq. (12.221) was derived in the p-regime of chiral perturbation 
theory in which mcL ^ 1 while mc^, c/L <^ p^. In particular, in eq. (l2.22p one 
cannot make Bg (and thus m) arbitrarily small, because one would otherwise enter 
the e-regime in which mcL ^ 1. 

Besides the constraint effective potential, Gockeler and Leutwyler have also de- 
rived analytic predictions for the first and second moment of the probability distri- 
bution p($), including the 2-loop level. The resulting expressions are 



(2.23) 
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For the cubic box considered here, the additional shape-coefficient is given by (32 = 
0.010608 p,]. 

Hasenfratz and Niedermayer have used the effective theory to derive the finite- 
size and finite-temperature effects of the staggered susceptibihty 



Xs 



1 + 2 



PsLl 



m) + 



PsLl 



1 



(2.24) 



from a 2-loop calculation in the e-regime of magnon chiral perturbation theory 
Similarly, the uniform susceptibility takes the form 



2p. 
3c2 



1 c ~ 1 
l + ^^/3i(/) + - 



3 PsLl 



3 \PsLl 



im - \ih{if - em 



(2.25) 

Here / = {Pc/LY^^ determines the shape of an approximately cubic space-time box 
of size L X L X /3, with /3c ~ L. The functions and i^i^l) are known 

shape-coefficients jsl, Isf. For an exactly cubical space-time volume (i.e. for / = 1) 
the result of eg. (12.241) agrees with eq. (l2.23p since 



(2.26) 



The factor 3 arises due to the three components of the staggered magnetization 
vector. 



3 Determination of the Low-Energy Parameters 

Ms, Ps, and c 

The susceptibilities Xs and Xu have been calculated numerically for the antiferro- 
magnetic spin | quantum Heisenberg model on the square lattice using the very 
efficient loop-cluster algorithm [10]. By comparing the Monte Carlo data with 
the effective theory predictions of eqs. fl2.24p and (I2.25p . the low-energy parame- 
ters Ms = 0.3074(4)/a^ = 0.186(4) J, and c = 1.68(l)Ja have been determined 
with high precision. At very low temperatures, one enters the cylindrical regime of 
space-time volumes with /3c ^ L in which the 5-expansion of chiral perturbation 
theory apphes. In this case, the staggered magnetization acts as a quantum mechan- 
ical rotor, again resulting in characteristic finite- volume effects |8|. By simulating 
cylindrical space-time volumes using the continuous-time variant of the loop-cluster 
algorithm, and again comparing with the corresponding predictions of the low-energy 
effective theory, the values of the low-energy parameters, previousl y o btained from 



the cubical space-time regime, have been verified independently HI]. Using the 



7 



value of the exact mass gap of the 2-d 0(3) model 2J], a result of Chakravarty, 
Halperin, and Nelson for the finite-temperature correlation length was extended 
by Hasenfratz and Niedermayer |25,] who obtained 



8 27rp, 



1 



+ 



(3.1) 



This expression is valid in space-time volumes with a slab geometry, i.e. for L ^ f3c. 
While the data in the cubical and cylindrical regimes determine the ratio Ps/c^ with 
higher precision than ps and c individually, Monte Carlo data for the very long cor- 
relation length in the slab regime determine ps = 0.1800(5) J very precisely 26|]. The 
combined analysis of all numerical data in cubical, cylindrical, and slab geometries 
resulted in c = 1.657(2) Ja and Ais = 0.30797(3)/a^ j26[. In a recent study using a 
zero-temperature valence-bond projector method, Sandvik and Evertz obtained the 
very accurate result ^As = 0.30743(l)/a^ Although the discrepancy between 
the two results for A4s is only in the permille range, it is statistically significant. In 
particular, it is important to clarify the discrepancy because our present very accu- 
rate study is sensitive to such small effects. For this purpose we have generated new 
data in the cubical regime for volumes ranging from L/a = 16 to 80. The largest 



volumes are substantially bigger than those of the original study [10|. In order to 
have an independent handle on the spinwave velocity c, in addition to the uniform 
susceptibility Xu, which is given by the temporal winding number Wt, we have also 
measured the spatial winding numbers Wi. The condition 



determines an exactly cubical space-time box with L 
til eq. fl3.2p is satisfied, we have determined c = L/ P = 



agreement with the result of [26 



(3.2) 

= l3c. By varying /3 un- 
1.6585(10) Ja in excellent 
By fitting the new data for Xs and Xu in the 
cubical regime to eqs. (12.241) and (I2.25p . we have obtained ps = 0.1808 (4) J and 
A4s = 0.30744(3)/a^. The value for ps is again in excellent agreement with the 
previous results obtained in the cubical, cylindrical, and slab regimes. The new 
result for Ais agrees with the one of the original study in the cubical regime l3 
and is consistent with the result of 27|, which is about two permille lower than the 



result of [26|]. We attribute this small but statistically significant discrepancy to an 



underestimation of the systematic errors of Xs in the cylindrical regime data of 26 
related to the termination of the Seeley expansion described in jsl. Until this issue 
is completely clarified, we discard the cylindrical regime data and instead include 
the result of The best estimate of the low-energy parameters obtained in cubic 
and slab geometries as well as at zero temperature is then given by 



Ms = 0.30743(l)/a^ = 0.1808(4)J, c = 1.6585(10) Ja. 



(3.3) 
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4 Improved Estimator for the Distribution of the 
Staggered Magnetization 



The loop-cluster algorithm j^, [l^, 11| is a very efficient numerical tool that al- 
lows us to perform high-accuracy numerical simulations of the quantum Heisen- 
berg model. The cluster algorithm connects the spin variables to closed loop-clusters, 
which are completely independent of one another. All spins in a given cluster are 
then flipped simultaneously with 50 percent probability. A given spin configuration 
containing N clusters is just one member of a sub-ensemble of 2^ equally probable 
configurations. An improved estimator substantially increases the statistics by ana- 
lytically averaging a given observable over all 2^ configurations in the sub-ensemble. 
For the quantum Heisenberg model, improved estimators have been constructed pre- 
viously for the stag gered and uniform susceptibilities Xs and Xu as well as for the 



energy density pj. The improved estimator for the distribution of the staggered 
magnetization to be constructed here is similar to the improved estimator for the 
topological charge distribution in the meron-cluster algorithm for the 2-d classical 



0(3) model 28|], which has been combined with a re- weighting technique 29 



In the loop-cluster algorithm, every cluster contributes additively to the total 3- 
component of the staggered magnetization. While it is straightforward to implement 
the improved estimator in continuous Euclidean time, it is most easily explained in 
the discrete-time variant of the loop-cluster algorithm 0, [ij. In that case, the 



cluster size \C\ (i.e. the number of lattice points in a given cluster) determines the 
3-component of the staggered magnetization of the cluster C, which is proportional 
to ±|C|. Under cluster flip the staggered magnetization of a cluster changes sign. 
The distribution of the staggered magnetization is recorded as a histogram which 
is built iteratively using one cluster after the other. The initial histogram ]9i(m) 
(with m G {— M, —M + 1, 0, M — 1, M}, where M is the number of space-time 
lattice points) is constructed from the first cluster as 

The two entries of the initial histogram correspond to the two possible orientations 
of the first cluster, each arising with probability |. In the z-th iteration step (with 
i G {2, 3, N}) a new histogram Pi{m) is built from the previous one as 

Pii^) = ^ [Pi-i{m + \Ci\) +pi^i{m - \Ci\)] . (4.2) 

After steps, all clusters have been incorporated and the final histogram is given 
by pAr(m). Examples of histograms pAr(m) obtained for two individual spin config- 
urations are shown in figure 1. The example in the left panel contains one cluster 
that is bigger than all the other clusters together. Hence, the region around m = 
is not sampled. In addition, there are two relatively large clusters that give rise to 
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m m 



Figure 1: Examples of histograms pN{m) obtained for two individual spin configura- 
tions on a 16^ lattice using the improved estimator. 

the multiple peaks in the distribution. In the example shown in the right panel, on 
the other hand, there are two clusters of similar size, such that the region around 
m = is also sampled. The final probability distribution of the 3-component of the 
staggered magnetization 

p(m) = (piv(m)), (4.3) 

is the average of the histograms pAr(m) for all configurations in the Markov chain 
generated by the cluster algorithm. By construction, this distribution is properly 
normalized as 

M 

P(^) = 1- (4-4) 

m=—M 

It should be noted that the numerical effort to build the improved estimator is 
proportional to the number of lattice points M and, in addition, proportional to the 
number of clusters. Since the number of clusters is proportional to the volume, the 
evaluation of the improved estimator requires a computational effort proportional to 
M^, and thus becomes rather time-consuming for large volumes. This is in contrast 
to the improved estimators for the susceptibilities Xs and Xu which only require a 
computational effort proportional to M. Of course, the improved estimator increases 
the statistics by a factor of 2^ which is exponential in the volume. Hence, investing 
a polynomial effort should still be justified. However, one should not forget that 
the 2^ configurations in a sub-ensemble are not statistically independent. As we 
will see, the improved estimator works very well and by far outperforms calculations 
done without it. 

The mean value of the 3-component of the staggered magnetization $3 corre- 
sponding to a given value of m is 

$3 = — • (4.5) 

^ 2M ^ ^ 
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The factor 2 in the denominator arises because we are deahng with quantum spins 
Hence, one can identify the probabihty distribution of the 3-component of the 
mean staggered magnetization as 

^($3) = p{m). (4.6) 



Due to the time-discretization chosen in [10|] (which has four Trotter steps) , all clus- 
ter sizes \C\ are multiples of 8. Consequently, the non-zero entries of the histogram 
p{m) correspond to values of m which are also multiples of 8. This implies that 



<i*s = ^. (4.7) 



Altogether, we thus obtain 



^, , p(m) M , ^ m , ^ 

p $3 = = — p m , $3 = , 4.8 

^ rf$3 4^^ ^' ^ 2M' ^ ^ 

with m restricted to multiples of 8. By construction, in the Euclidean time contin- 
uum limit the resulting probability distribution is normalized as 

rf$3 p(<f 3) = 1. (4.9) 



We have simulated the spin | quantum Heisenberg model on a square lattice 
with L/a between 8 and 24 at inverse temperatures (3 = L/c using the loop-cluster 
algorithm in its discrete-time variant j^, [lo[. We have worked at a sufficiently small 
lattice spacing in Euclidean time, such that the systematic discretization error is 
negligible compared to the statistical errors. The probability distribution ^($3) of 
the 3-component of the mean staggered magnetization $3 has been obtained using 
the improved estimator described above. A typical distribution is shown in figure 
2. As we will see below, the information about the vicinity of the minimum of 
the constraint effective potential «($) is contained in the region of $3 where ^($3) 
changes rapidly. Figures 2 and 3 compare Monte Carlo data obtained with and 
without the improved estimator, investing the same amount of computer time in 
both cases. The error reduction of the improved estimator is very substantial. By 



employing the re- weighting technique used in [28|, 129|, one could concentrate the 
statistics in the relevant region in which p($3) changes rapidly. This should further 
improve the efficiency of our numerical method. 

Due to the 0(3) symmetry, for a fixed magnitude $ of the mean staggered 
magnetization vector $, its 3-component has a fiat distribution given by 

M^^) = ^ 0*(*3). (4.10) 

Here 0<j)($3) is a step function which is equal to 1 for $3 G [—$,$] and zero 
otherwise. The probability distribution of the 3-component ^($3) and the one of 
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Figure 2: Probability distribution p{^z) of the 3-component of the staggered magne- 
tization $3 on a 16^ lattice obtained with the improved estimator. The error bars of 
the distribution are on the order of the line width in this figure. 




Figure 3: The probability distribution ^($3) of the 3-component of the staggered 
magnetization $3 on a 16^ lattice obtained without the improved estimator. Unlike 
in figure 2, which contains more than 4x10^ bins, here the data have been compressed 
into 10^ bins, which substantially reduces their otherwise even larger variance. 
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the magnitude are related by 



p($3)=47r/ d$ $2 p$($3) = 27r / rf<l> $ p(<l>) 0$(<l>3), (4.11) 
Jo Jo 



such that 



/■oo 

27r / rf<l> $ J9(<1)) (5(<1>3 - $) = -27r<l)3 j9($3). (4.12) 
Jo 



Hence, given the Monte Carlo data for ^($3), we can extract the probabihty distri- 
bution of the magnitude of the staggered magnetization as 



47r$' p($) = -2$ ^^4^. (4.13) 
This indeed ensures the correct normahzation of eq. fl2.10l) because 

47r [ d<l> $2 p($) = -2 [ rf$ $ = 2 [ d<^ p{<^) = [ rf$3 ^($3) = 1. 

Jo Jo Jq 

(4.14) 

Using eq.f l4.13p . we have determined the probabihty distributions p($) from ^($3) 
obtained using the improved estimator. Some results for 47r^$^p($) are shown in 
figure 4. As the volume increases the mean value of $ decreases and the width of 
the distribution p($) becomes narrower. It should be noted that the distribution is 
not symmetric around its maximum. Injiie infinite volume limit, the distribution 
turns into a 5-function centered at $ = J^s = 0.30743(1). 

We also like to compute the first and second moment ($) and (($ — ($))^) of 
the distribution p($). Of course, this is trivial once we have computed p($) us- 
ing the improved estimator which requires a computational effort proportional to 
M^. If one only wanted to compute the two moments but not p($) itself, one may 
wonder how to achieve this in the most efficient manner. As we noted before, the 
combination (($ — ($))^) + ($)^ = ($^) is proportional to the staggered susceptibil- 
ity which can be obtained using an improved estimator requiring a computational 
effort proportional to M. How can one determine ($) itself? Since $ is the magni- 
tude of the mean staggered magnetization vector of which only the 3-component 
is easily accessible in a quantum Monte Carlo simulation, this seems not entirely 
straightforward. However, using eq. fl4.13p one obtains 

($) = 47r^ p{^) = -2 c;$$2^^^1 = 4^ rf<l> $ p($) 

/oo 
t/$3 |$3|p($3)=2(|<l>3|). (4.15) 
-00 
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Figure 4: Probability distributions 47r$^p($) of the magnitude of the staggered mag- 
netization $ for L = 8a, 12a, 16a, 20a, and 24a. The vertical line at $ = Aig = 
0.30743(1) represents the 6-function distribution of the infinite system. 

Hence, by just measuring |$3| in a standard unimproved manner, one can determine 
($) with a computational effort proportional to M. The crucial question is how the 
statistical error achieved in this way compares with the one obtained by extracting 
($) from (whose construction requires a computational effort proportional to 
M^). In order to investigate which of the two methods is more efficient, we have first 
determined ($) from p($) using the improved estimator described above. Then we 
have invested the same amount of computer time in an unimproved measurement 
of 2(|$3|). For L/a = 20 and f3 = L/c the statistical error of ($) obtained in this 
manner is a factor of 1.5 larger than when one uses the improved estimator. Hence, 
despite its computational effort proportional to M^, thanks to the average over the 
large number 2^ of configurations in the sub-ensemble, the improved estimator is 
slightly better than the unimproved method, even if one is only interested in the 
first moment ($) and not in the distribution p($) itself. 

In order to compare our Monte Carlo data with the expectation value |($)(i?s)| 
of eq. fl2.22p . we also like to switch on a staggered magnetic field Bs = (0, 0, Bg). The 
corresponding probability distribution 

P($3, Bg) = -^fi^s) exp [^^^^) ■ (4.16) 

can be sampled using the loop-cluster algorithm with an additional Metropolis 
accept-reject step applied to each cluster flip, which takes into account the contri- 
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bution exp($3i?5L^/5/a^) to the Boltzmann weight. Using an ordinary unimproved 
estimator, one then simply measures = |($3)|. 

5 Comparison of Monte Carlo Simulations and 
Effective Theory Predictions 

We have determined the first and second moment ($) and (($ — ($))^) of the dis- 
tributions which are compared with the effective field theory predictions of 
eq. fl2.23p in table 1. The errors of the theoretical predictions result from the uncer- 
tainties in the low-energy parameters of eq. fl3.3p . For the first moment the agreement 
is very good for L/a> 10. The absolute value of the second moment is very small 
and its statistical error is relatively large. Still, there are systematic discrepan- 
cies between the Monte Carlo data and the (9(1/L^) effective theory predictions of 
eq. fl2.23p . This discrepancy is well accounted for by additional 0{1/ Li^) corrections. 
Such corrections involve new low-energy parameters multiplying higher-order terms 
in the effective action. Their evaluation would require a 3-loop calculation which 
has not been worked out in the effective theory. Parameterizing the 3-loop terms 
with unknown coefficients ai and a2, i-e. 

((4-(*>r>-f^A + «.(^)%o(ij). (5.1) 

one obtains a good fit to the Monte Carlo data for ai = —0.0017(2) and 0^2 = 
—0.00042(2). One may conclude that precise calculations of the two moments allow 
the determination of some combination of sub-leading low-energy parameters. 

Using p($) = A/'exp(— L^m($)), the probability distributions of figure 3 are 
readily converted into the corresponding constraint effective potentials m($) shown 
in figure 5. With increasing volume the constraint effective potential approaches the 
convex shape of thejnfinite volume effective potential. Using the rescaled variable 
ip = {psL/c){^ — A4s)/-M.s, one can also consider the extensive quantity Ui^ip) = 
L^m($) which is shown in figure 6. Expanding f/($) = Uo{ip) + (c/psL) Uiijp) + 
0(1/L^), we have used the Monte Carlo data for L/a between 12 and 24 to determine 
Uq{i{j) and Ui{ip). Some values of the function Uo{ip) extracted from the numerical 
data are compared with the analytic result of eq. 02.161) in figure 7. It should be 
pointed out that the observed agreement does not rely on any adjustable parameters. 
Even the normalization constant A/" of eq.f l2.12l) . which fixes an additive constant in 
the constraint effective potential, is predicted by the effective theory. As quantified 
in table 2, in the interval G [0.10, 0.35], i.e. around the minimum of the constraint 
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T In 

Li 1 a 




\^/ theory 


- ;<P;j ;mc 


- (<I>)j )thcory 


8 


0.38841(17) 


0.3913(2) 


7.1(3) xlO ^ 


1.318(6) xlO 


10 


0.37258(14) 


0.3738(2) 


5.1(3) xlO ^ 


8.44(4) xlO * 


12 


0.36186(13) 


0.3624(1) 


3.6(3) xlO * 


5.86(3) xlO * 


14 


0.35401(12) 


0.3543(1) 


2.9(2) xlO * 


4.30(2) xlO ^ 


16 


0.34810(8) 


0.3483(1) 


2.5(2) xlO ^ 


3.30(2) xlO * 


18 


0.34339(13) 


0.3437(1) 


2.3(2) xlO 


2.60(1) xlO 


20 


0.33990(9) 


0.3400(1) 


1.6(2) xlO-^ 


2.11(1) xlO"^ 


22 


0.33683(15) 


0.3369(1) 


1.6(2) xlO-^ 


1.74(1) xlO-^ 


24 


0.33462(19) 


0.3344(1) 


1.0(3) xlO-^ 


1.46(1) xlO-^ 



Table 1: Comparison of Monte Carlo data (MC) for the first and second moment 
($) and (($ — ($))^) of p{^) with the effective theory predictions of eg. h2.2^) at 
the 2-loop level. The errors of the theory predictions are due to small uncertainties 
in the values of the low-energy parameters of eg. ^3.3^) . The discrepancy between the 
Monte Carlo data and the effective field theory results is due to a 3-loop correction 
that was neglected in the theoretical predictions. 




0.44 



Figure 5: Constraint effective potentials «($) as functions of the magnitude of the 
staggered magnetization $ for L = 8a, 12a, 16a, 20a, and 24a. The constraint effec- 
tive potential approaches a convex effective potential in the infinite volume limit. 
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Figure 6: The extensive quantity U{^) as a function of the rescaled variable ip = 
(psL/c)($ — J^s)/-M.s for L = 8a, 12a, 16a, 20a, and 24a, compared to the analytic 
infinite volume result Uolip). 




Figure 7: The analytic result for the universal function Uo{ip) compared to numerical 
values obtained from a fit of the Monte Carlo data for f/($) to eg. I\2.14\ l- 
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f^o(^) theory 


0.10 


- 3.782(62) 


- 3.818 


0.15 


- 4.271(39) 


- 4.302 


0.20 


- 4.578(22) 


- 4.589 


0.25 


- 4.646(18) 


- 4.653 


0.30 


- 4.458(20) 


- 4.462 


0.35 


- 3.983(36) 


- 3.986 



Table 2: Comparison of Monte Carlo data (MC) for the universal function Uq{iI)) 
with the effective theory prediction of eqs A2.1b\) . 



L/a 


Bs/J 


\mBs)\ 


24 


0.0300 


0.35269(1) 


24 


0.0330 


0.35472(1) 


32 


0.0175 


0.34302(1) 


32 


0.0185 


0.34396(1) 


40 


0.0100 


0.33487(1) 


40 


0.0110 


0.33615(1) 


48 


0.0075 


0.33142(1) 


48 


0.0080 


0.33219(1) 



Table 3: Monte Carlo data for |($)(i?s)| which are used in the determination ofko. 

effective potential, the theoretical values and the simulation data for Uo{ip) agree 
within error bars. 

In order to determine the low-energy parameter fco, we have considered the expec- 
tation value |($)(-Bs)| as a function of the external staggered magnetic field given in 
eq. fl2.22p . Since these results do not require the probability distribution p($), they 
could be obtained on larger volumes up to L/a = 48. The results summarized in 
table 3 give a good fit for 

fcg = -0.0037(3). (5.2) 

Using the theoretical prediction for Uo{ip), figure 8 compares values of the func- 
tion Ui{ip) determined from the Monte Carlo data with the analytic result of eq.f l2.18p . 
for the two values of ko at the edge of the corresponding error band. In the interval 
ip G [0.125, 0.3] around the minimum of the constraint effective potential, the Monte 
Carlo data are consistent with the theoretical predictions of Gockeler and Leutwyler. 
We attribute the deviations outside this interval to effects of higher order. Indeed, 
as one sees in figure 6, the finite volume effects become larger when ijj moves away 
from the minimum of the constraint effective potential. 
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Figure 8: The analytic result for the function Ui{il)) compared to numerical val- 
ues obtained from a fit of the Monte Carlo data for U{(^) to eg. 1^2.14 ), using the 
theoretical prediction of Uq{iP) as an input. 



6 Conclusions 



We have constructed an improved estimator for the probabihty distribution of the 
staggered magnetization in the quantum Heisenberg modeL Using the improved 
estimator in a loop-cluster algorithm simulation, we have determined the first and 
second moment ($) and (($ — ($))^) of the distribution of the magnitude $ of 
the mean staggered magnetization vector $, as well as the constraint effective poten- 
tial u{^) (obtained from = A/'exp[— for different space-time volumes. 
The Monte Carlo data are in excellent quantitative agreement with analytic predic- 
tions which Gockeler and Leutwyler derived from a systematic low-energy effective 
field theory. This demonstrates that the magnon effective theory indeed provides 
correct predictions, order by order in a systematic low-energy expansion. Thanks to 
the very efficient loop-cluster algorithm, the Heisenberg model is an excellent testing 
ground for the effective field theory method. Its quantitative success provides en- 
couragement to also test the recently constructed systematic effective field theories 
for holes 3^, 31] and electrons 32] doped into an antiferromagnet against numerical 
simulations. First results confirming the effective theory have already been obtained 



for the t-J model on the honeycomb lattice ^2|. The results obtained in the study 
presented here should also be encouraging for lattice QCD simulations, were the 
numerical problem is much harder. Eventually, one may expect agreement between 
lattice QCD and chiral perturbation theory at the same level of accuracy as achieved 
in the condensed matter problem investigated in this paper. 
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